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SUMMARY 


Lagrange's method of undetermined multipliers is applied to the 
problem of properly combining lift loadings for the least drag at a given 
lift on supersonic wings. The method shows the Interference drag between 
the optimum loading and any loading at the same lift coefficient to be 
constant. This is an integral form of the criterion established by 
Robert T. Jones for optimum loadings. 

The best combination of fo\ar loadings on a delta wing with subsonic 
leading edges is calculated for several Mach numbers as a numerical exam- 
ple. The loadings considered have finite pressures everywhere on the 
plan form. At each Mach number the optimum combination of the four non- 
singular loadings has about the same drag coefficient as a flat plate 
with leading-edge thrust. 


INTRODUCTION 


The problem of minimizing the supersonic drag for a given lift on a 
fixed plan form has been approached in different ways. ■ Jones, in refer- 
ences 1 and 2, makes use- of reverse-flow theorems to derive several simple 
properties of the optimum load distribution and to present as well the 
optimum distribution for elliptic, plan forms. Graham, in reference 5, 
shows how the proper use of orthogonal loadings can reduce the drag at 
fixed lift. Orthogonal loadings .are loadings of zero interference drag. 
The interference drag between two loadings is the total drag of each in 
the downwash field of the other. In reference 4, theorems concerning 
orthogonality and reverse flow are developed, whereas in references 5 
and 6 mmierical examples of drag reduction by use of orthogonal loadings 
are given. For delta wings with conical camber the optimum shapes are 
derived by Ritz's method in reference 7 . 

In this paper Lagrange's method of undetermined. multipliers is 
applied to the problem of properly combining loadings for the least drag 
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at a given lift. By use of this method a simply expressed property of 
the optimum loading is found which is an integral form of a property 
established by Jones in reference 1. Jones' property of the optimum 
loading is that the downwash on the plan form is constant in the combined 
forward- and reverse-flow fields. The best combination of four types 
of nonsingular loading on a delta wing is calculated as a numerical 
example of the use of the method. 


A 

b 

c 


SYMBOLS 

loading strength parameter 
span 

local chord 
drag coefficient 


*^D,ii ^D,i drag coefficient, of ith loading 

Cj) 4 j drag coefficient of interference between ith and jth 

^ component loadings 

Cl lift coefficient 

^ lift coefficient of 1th loading 


Cj(x) lift coefficient on spanwlse line 

Cj (y) lift coefficient on local chord 

Cp lifting pressvire coefficient 

M Mach nvanber 

m tangent of semiapex angle 

N nvanber of loadings 


n = pm 


R 


functions of 0 and n (see appendix) 


S 


wing area 
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s,t 

x,y,z 

a 

3 = Vm^ - 1 

€ 

0 

mx 

A 

T 

1 - 

Subscripts: 

i,J 

M 

0 

X 


integers 

Cartesian coordinates of lifting surface (see fig. 2) 
local angle of attack of lifting surface 


small positive number 


Lagrange's multiplier 
plan form 

root chord of arrow wing 

ith, jth loading component 
minimum among all loadings 
minimum among N loadings 
arbitrary loading 


ANALYSIS 

Theory 


Consider a superposition of N loadings of the form 


- ■^iSa ^3^A + • • • + ^N^p,N 

where A is the strength parameter and Cp is the resultant lifting 

pressure coefficient at a point on the plan form. The corresponding 
local angle of attack may be written as 


a = + Agttg + AjO^ + . 


■ • Vn 


( 2 ) 
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The local drag coefficient Cpa is a quadratic in A which may be 
integrated over the plan form t to give the drag coefficient of the 
wing. Thrust-loaded singularities at the leading edge are therefore 
excluded from the drag. This exclusion is merely for convenience and 
is not necessary. The formula for the drag coefficient is 

- 1 / ^ ^ (5) 

i=l J=i . 

The average lifting pressure coefficient on the plan form is the 
lift coefficient, which is 

r ^ 

Cl - I / Cp ds . X Cl,iAi P*) 

i=l 


The problem is to find the set of A's which yields the minimum 
value of Cp subject to the condition that Cp is constant. Because 
of the quadratic nature of Cp and the linear form of Cp, Lagrange's 
method of ixndetermined multipliers is particularly suitable for the, 
solution as it leads to a set of linear algebraic eqiiations. 

As shown in reference 8, a function of the A coefficients 
F = Cp + ACp is formed, where 'K is Lagrange's multiplier. The minimum 

value of F as determined by the N linear algebraic equations = 0 

oAp 

plus condition (4) is Lagrange’s solution. A schematic drawing of the 
coefficient matrix of the N + 1 linear equations is shown in figure 1. 


The equations may be written more simply if first the interference 
drag between the optimum loading and the ith component of the loading 
is computed. From equations (l) and (2), the following expressions may 
be written: 


^p,0°'i = ^3^P,5°1‘‘' -^N^p,N^i 


(5) 


Adding equations ( 5 ) and integrating over the plan form gives 

^D,Oi = lX(S,0°^i + + ^ 2%, 21 + ^5'^D,5i 


• • • + 


^^i^D,ii +•••■*■ ^N^D,iN 


(6) 
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This expression for Cp Qj^ is a part of the left-hand side of the ith 
equation of the linear set which is now written as 

^D,Oi ° ('7') 

A simple property of the optimum load distribution may now be derived. 
First written in terms of" 


^D,0 = •^l(^D,01 " ApCD^ii) + ^2{^’D,02 " ■^1^D,12 ~ ■'^2^D,22) + • • • + 

^n[^D,ON ' ■'^1^D,1N ■ -^2^,2N " '^5*^0, 5N ’ • • • ' 

The sum of the terms not containing is -Cj)^0* Substituting for 

Cjj Qj^ from eqiiation (j) gives 


^D,0 - "‘^D,0 


'(^ 1 ^ 1,1 + ^ 2 ^ L ,2 + ^ 3 ^ L ,3 + • • • + 




(9) 


The term in parentheses is Cj^; therefore A may be written as 


D,0 


A = -2 


( 10 ) 


Substituting equation (lO) into equation (j) gives 


n _ o p 

CD,0i - 2 Cl,1 


( 11 ) 


Since equation (ll) holds for any number of loadings, let the number 
of components increase without limit to include all possible loadings. 

For an arbitrary loading X and the absolute minimxmi M, equation (ll) 
may be written as 

C ' 

^,MX = 2 Cl^X (12) 

The meaning of equation (12) may be simply expressed as follows; 

The interference drag between the optimum loading and any loading at the 
same lift coefficient is constant. If the reversibility theorem is applied, 
equation (12) is an integral equivalent of a condition established by 
Jones in reference 1. Jones' .condition states that for the optimum loading 
the downwash on the plan form is constant in the combined forward- and 
reverse-flow fields. Barred variables will represent the reverse flow 
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which has the same lift loading on the plan form hut, in general, a 
different sxirface shape. Then, by reversibility, 

By definition, Cp^MX is 

^D,MX "sf °M^p,x)‘^^ 


Therefore^ equation (12) may be written as 







Since Cp^x 


is arbitrary, + % must be constant. ■ Hence, 


0^ + 5^ 



( 14 ) 


( 15 ) 


This is the condition derived by Jones in reference 1. Equation (12) 
is then an equivalent integral form of this condition. 

Eqiiation (12) shows the orthogonality of the optimum loading to, 
and only to, zero lift loadings. This result, which was stated by Graham 
in reference 5, is seen to be a special case of a more general inter- 
ference drag property given by equation (12). 


Comparison With the Method of Orthogonal Loadings 

If two loadings are to be combined, it may be shown that Graham's 
method of orthogonal loadings (ref. 3) and the present method are 
equivalent. If the resultant combination of two loadings is combined 
by the method of reference 3 with a third loading, the lift ratio of 
the first two loadings is unchanged in the best combination of the three. 
If n > 2 loadings are successively combined in the manner of refer- 
ence 3, the first n - 1 loadings are not allowed to adjust their 
relative strengths upon addition of the nth. In the present Lagrangian 
method every loading has equal freedom to adjust. For this reason, the 
Lagrangian method should be more rapidly convergent. 



NACA TN 5535 


7 


NUMERICAL EXAMPLE 


Tucker in reference 9 presents formulas for the surface coordinates 
of delta and arrow wings which support four types of pressure distri- 
bution. The formulas axe given for subsonic leading edges and supersonic 
trailing edges. In the notation of this paper (fig. 2) a combination of 
the four loadings may be written. 

Cp = Ai + A2X + Aj + A4 ^ (16) 

Formulas for the QD,iJ quantities may be derived from equation (I6) and 
the surface formulas given in reference 9> Ly integrations over the plan 
form. Details are given in the appendix. 

The optimum drag results and strength parameters are presented in 
figures 5 and 4. In figure 5 the corresponding drag values for a flat 
delta wing with and without leading-edge thrust (ref. 10) are given for 
comparison. The drag values for the four component loadings taken alone 
are also shown. In addition^ the drag of the conically cambered optimvim 
delta wing (ref. j) aind Jones' absolute minimum for narrow wings (ref. 1) 
are plotted. 

Noteworthy in figure 5 is the closeness with which all the optimum 
drags agree with each other and with the drag of a flat delta wing which 
has a thrust-loaded leading edge. The close approach of the present 
optimum of four loadings to Jones' absolute minimum for narrow wings is 
also evident . 

In figure 5 the chordwise and spanwise loadings of the optimum 
combination are shown for the extreme ^ch numbers calculated. Compar- 
ison with the elliptical loading shows good agreement spanwise and poor 
agreement chordwise. 

The data indicate that the relatively low drag of the flat delta 
wing with leading-edge thrust can be equalled by properly combining a 
few loadings having finite pressures everywhere on the plan form. 

A plausible specvilation suggested by the data is that it is possible 
to come very close to the minimum drag on a delta wing with but a few 
steps in a series approximation. Perhaps, too, a restricted minimum, 
such as the one for conical camber, gives a close approximation to the 
absolute minimum drag if the restriction is not too unnatural. 
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CONCLUDING REMARKS 


Lagrange's method of undetermined multipliers is applied to the 
problem of properly combining lift loadings for the least drag at a 
given lift on supersonic wings. 

The method shows the interference drag between the optimton loading 
and any loading at the same lift coefficient to be constant. This is an 
integral form of the criterion established by Robert T. Jones for optimum 
loadings . 

The best combination of four loadings on a delta wing with subsonic 
leading edges is calculated for several Mach numbers as a numerical 
example. The loadings considered have finite pressures everywhere on 
the plan form. At each Mach number the optimijm combination of these 
four nonsingular loadings has very nearly the same drag coefficient as 
a flat plate with leading-edge thrust. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 2T , 1955* 
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APPENDIX 

DETAILS OF NUMERICAL EXAMPLE 


Inasmuch as the pressure coefficient and corresponding a are 
given by 

|y I 

Cp = Ai + A2X + A3 -1^ + A4 — 

O' = ApO^ + Apttp + A^o^ + 

then the local drag coefficient may be written as 

Cptt = + AiA2(xai + 02 ) + AiA3^-i^ai + + 

^2^ (^“2) ^2^5 (^^“2 ^“'5^ 

{^4 


+ xoi^ + k-^ 

4^ 


^ ^3A4(^cc5 


(Al) 


^2^5 + 


(A2) . 


The required functions are the averages over the plan form 

(fig. 2) of the quantities in parentheses in equation (A2). Rather than 
tti itself, reference 9 gives the surface ordinate which is the 

chordwise integrated value of a^; 


Zi = - / ai dx 


The values given for are 


x2 

^2 = -sr ^2 

x2 

^5 = - 1^5 

x5 

^4 = - ^4 




(A3) 


(A4) 
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The values of R^ are functions of 0 = tabulated in reference 9 

mx 

for different values of n. The equations for are 


R, = r^2^(l~-~TF^ - 2 cosh"^ 1-^1 + - n^(l + 0) cosh”^ 

1 4it |_ I n0 1 


1 + n^0 


n(l + 0) 


- n^(l - 0) cosh"^ 


1 - n^0 


n(l - 0) 


(A5a) 


R2 = 




)Jl - rf 


\jl - n202 _ 202 cosh-1 + 


n^(l ~ ^ ^ + 0 + ©2 I cosh"^ 


2 


1 + n^0 


n(l + 0) 


0 02] eo=-h-^ 

H 

1 

CD 


L 2 J 

CD 

1 

r — 1 



(A5b) 


R^ = -r^ 


1_ 

4n 


I ^1 - n^02 - + 502 _ I n202j cosh"^ 


_ 1 ^ 

n0 


(1 + ef + 2(1 - n2)(0 + 02) 

1 + n20 


(1 !“o 

n(l + 0) 


2\Jl - n^ 



(1 - 0)2 - 2(1 - n 2)(0 - 02 ) ! 

1 - n2e J 


2/1 - n2 

n(l - 0) 



(A5c) 
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(l - n2) 


(l - . 12 - lOn^ 2n2\l^ 2^2 ^n2 ^-1 

^ 7 oV + -~^/ o\ ^ U - n'^^e - d 6 cosh 

3(1 - n2) 3n2(l - n^j 

2n“(,2 ^ ,5) ^ 2^{e . e 3 ) . 


n0 


W 


f (1 . e3)_ 


£f 

ln(l + 0) 


cosh-1 l _ t n^e ^ 2^ 


(1 - n2) 


3/2 


|6 - 9a|,+,. 2n^(e2 _ e3) . 


(0 - ©5] 

_ £(i . 05) 

.-1 

cosh 

i-„^e 1 1 


6 ' . 

J 

CD 

1 


(A5d) 


For terms in equation (A2) of the type (y/m)^OQ_, a spanwise integration 
of gives the following average on the plan form: 




Rid) ,, ,s+t+l e®R^(0) _ 

. io ' ' 


For terms of the type xaj^ an additional integration by parts in 
the X direction is required to maintain the Rj|^ functions intact under 
the integral signs. The result for this case is 


i / xai' dS = 


— 2__, 1 
1 - m 


%d) ,t+2 pi R.(0) 

- d - ^i) / i rr^d0 + 


0 (1 - ^i0)l■^^ 


(1 - Rj(0) 


iO r 

+ 2 Jo 


d - n0) 


t+2 


d0 


(A7) 


In formulas (a 6) and (AT) the value of t for each 1 is as follows: 


1 

2 

3 

4 


1 

2 

2 

3 
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By applying formulas (a6) and (A 7 ) to the integration of (A2), the 
followi^ equations for Cp are derived: 


^%,11 = ^ 1^1(1) - ^(1 - ki) 



Ri(e) 

(1 - pe)5 


de 


(A 8 a) 


"Cd ,12 = ^ 


2 ^1 


Ri(e) 


de - 


2(1 - hr 


Ri.(e) 


0 (1 -,h0) 


de - 2(1 - 


'0 (1 - ^,0.)^ 
2 pi Rp( 0 ) 


0 (1 - ne)^ 


de. 


(A8b) 


“•^,15 = 3 ~ (T~u) - 2(1 - h)^ r de - 

^ Jo (1 - ne)^ 


2(1 - h)‘ 


2 . %(®^ 

(1 - ne)^ 


'f 

J 0 


de 


(a8c) 


“>Cd,i 4 = ^,) [gl(l) + %(1^ - 2(1 - h)^ 

2(1 - 

■ J n 


l^i e^R^Ce). 

0 (1 - ^ 6)5 


de - 


0 , ( 1 . - ne)^ 


(A 8 d) 


2mC-Q^22 


RpCl) 

1 - h 


-4(i-u) 5 r" 
J 0 


R 2 (e) 

(1 - h0)^ 


de + (1 



R 2 ( 0 ) 

— £ de 

(1 - h0)^ 


' (A 8 e ) 


NACA TN 5553 


13 


1 

2(1 - n) 


|r2(1) + R5(1]) - 2(1 - _ 

-^0 (1 _ ne )5 


Jo (l - 110)^ 


H^(e) 


'Q (1 - 110) 




R4(6) 


0 (1 - n0)° 

1 R4(6) 


2(1 -n)M - . -t . / ae-f^d-n)^ — ae 

Jo (1 _ ^, 0)6 5 Jo (I . 




1 - \1 


0 (1 - H0)5 


>>3** ' '*■>(1^ - 2(1 - ^‘)‘ 


'1 02R^(e) 

> (1 - H0)' 


2(1 - 44)^ f 

J n 


0R^(0) 


0 (1 - ki6) 


^ ^0 (1 _ ^0)7 


1 9%i^(0) 


Ik 
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The required 
form which yield 


"L,i 


functions are simple integrals over the plan 




(A9) 


The integrals in equations (a 8) were, in general, evaluated numer- 
ically. However, several of the integrands in equations (a 8) have the 

RiCe) R.Ce) 

form — and — r-. These functions have an infinite 

(1 - H6) (l - 

discontinuity at 6=0. For such a discontinuity, numerical methods 
break down. Near zero the following approximation is integrated 
analytically: 


Si(e) « Ei(e) .. i oosh-l ^ i cosh-1 i 
Rj(9) . Rj(,) . i cosh-1 i ^ i i 


W0< 0'S €« 1 


(AlO) 


The integrals for the region 0 = 6 = e can then be approximated: 


Rl(0) 


'o (1 - lae)* 


de 


f(s)[Ei(c) ^ ^ cosh-1 i] - 


de « f(e)lR,(e) - 5-^ cosh"^ J 
- ' Art nej krt 


'0 (1 - iiQ)^ 


^ (All) 


where 


f(e) = 


de 


'0 (1 - p 6 ) 


= e 


_2 


^ t e t(t + 1) pe t(t -f- 1) (t -1- 2) US' 
2: 3 3: A 


+ • . • 


(A12) 
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and 


:(.) =/ 

'^0 


e cosh-1 ^ 

(1 - n0)* 


d0 


(A15) 


The integral in equation (A13) may be evaluated by expanding the 
denominator by the binomial theorem and writing l(e) as an infinite 
series 


l(e) = a^iQ + a^i^ + a^i 2 + a^i^ + 


(A14) 


where 

ao = 1 
t 

ai = 3 - h. 


a.2 ” 


_ t(t + 1) 2 


21 


lo-^ cosh-1 i 
' / 


d0 


0 cosh“l — d0 
n0 


t(t -h l)(t 2) 
a-z - 

^31 

t(t 1) . . . (t -H S - 1) s 
as - ^ ^ h 


-r 

'-^0 

-f 

^0 


r (A 15 ) 


0 ^ cosh 


"1 d0 
n0 


0® cosh“l d0 
n0 


The ig integrals of equation (AI 5 ) are evaluated by use of the 
relation 


cosh-1 Ji. (10 = : cosh-1 ^ __1 — r Q 

J n 0 s-hl n 0 s-flj 1/ 2^ 

Vl - n 0 


d0 (a16) 
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Figure A.- Optimtun values. of A parameters on a delta plan form. 








